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ABSTRACT 

We analyze the critical behaviour of the 3-d, 3-state Potts model in the presence of 
an external ordering field. From a finite size scaling analysis on lattices of size up to 
70'^ we determine the critical endpoint of the line of first order phase transitions as 
{/3^,hc) = (0.54938(2), 0.000775(10)). We determine the relevant temperature like 
and symmetry breaking directions at this second order critical point and explicitly 
verify that it is in the universality class of the 3-d Ising model. 
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1 Introduction 



The temperature driven first order phase transition of the three-dimensional, three- 
state Potts model has been analyzed in great detail in the absence of a symmetry 
breaking external field [||, 0]. This transition remains first order in the presence of a 
non-vanishing external field and for a critical field strength, he-, it ends in a second 
order critical point. Although it is expected that this critical point belongs to the 
universality class of the 3-d Ising model, the universal behaviour in its vicinity has 
not been analyzed in detail so far. A first estimate for the location of the critical 
endpoint is given in |]^. 

Our interest in properties of the 3-d, 3-state Potts model in the vicinity of the 
second order critical point is motivated by its importance for the analysis of lines of 
first order transitions in lattice gauge field models. The universal behaviour in the 
vicinity of the second order endpoint of the line of first order transitions has recently 
been investigated in the U(l)-Higgs [Q and SU(2)-Higgs |]^ models. Lines of first 
order transitions with critical endpoints do, however, also occur in QCD with light 
as well as heavy quarks. In the heavy quark mass limit of finite temperature QCD a 
line of first order phase transitions (deconfinement transition) occurs which is closely 
related to the phase transition in the 3-d, 3-state Potts model in an external field 
P, 0. Like in the Potts model the deconfinement transition is first order in the 
limit of infinitely heavy quarks and ends in a second order transition at some finite 
value of the quark mass. This critical point is expected to belong to the universality 
class of the 3-d Ising model [Q. Another line of first order transitions occurs in the 
case of QCD with three light quark flavours (chiral symmetry restoration). This 
line also ends in a second order endpoint at some critical value of the quark mass. 
In order to analyze the critical behaviour at these endpoints it will be important to 
disentangle the relevant energy- and ordering field like directions and identify the 
Ising-like observables at the critical points. We will address this problem here first 
in the simpler case of the Potts model and will explore methods used for the analysis 
of the liquid-gas phase transition as well as lattice gauge models 0, 

In this letter we present an accurate determination of the critical couplings at the 
endpoint of the line of first order phase transitions in the 3-d, 3-state Potts model 
with an external ordering field. Moreover, we will determine the relevant energy-like 
and ordering field like directions at this endpoint as well as the related operators 
from which critical exponents and other universal constants (Binder cumulants) 
can be extracted. In the next section we fix our notation for the Potts model and 
introduce the new couplings and operators which control the critical behaviour at 
the endpoint. In section 3 we present our numerical results for the location of the 
critical endpoint. A more detailed discussion of the universal properties at this 
endpoint is given in Section 4. Finally we present our conclusions in Section 5. 
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2 The Model and Simulation Parameters 



The three-state Potts model is described in terms of spin variables (Tj G {1, 2, 3}, 
which are located at sites i of a cubic lattice of size V = L"^. The Hamiltonian of 
the model is given by, 

H = -I3E -hM , (2.1) 

where E and M denote the energy and magnetization, 



E = Y,6{a,,a,) , M = Y,6{a,,ag) . (2.2) 

(id) * 



Here the first sum runs over all nearest neighbour pairs of sites i and j. A non- 
vanishing field h > favours magnetization in the direction of the ghost spin (jg. On 
a finite lattice of size the partition function of the model is then given by 

Z(/?,/i,L) = ^e-^ . (2.3) 

For vanishing external field the model is known to have a first order phase transition 
for j3c{h = 0) = 0.550565(10) 0. In the presence of a non- vanishing external field 
h this first order transition weakens and ends in a second order critical endpoint, 
which is expected to belong to the universality class of the 3-d Ising model. A first 
estimate of the critical endpoint (/3c, has been given in Ref. 0. We will give here 
a more precise determination of (/3c, ^c) and analyze the universal critical behaviour 
at this point. 

At {(3ci he) the original operators for the energy and magnetization, E and M, 
loose their meaning as operators being conjugate to the temperature-like and sym- 
metry breaking couplings. One rather has to determine the new relevant directions 
at the critical endpoint, which take over the role of temperature-like and symmetry 
breaking directions and allow the determination of the two relevant critical expo- 
nents at the second order endpoint. This also fixes the new order parameter and 
energy-like observables as mixed operators in terms of the original variables E and 
M. Following the discussion of the corresponding problem for the liquid-gas transi- 
tion we introduce new operators 

M = M + sE , E = E + rM , (2.4) 
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as superpositions of the original variables E and M. The Hamiltonian of the Potts 
model can then be rewritten in terms of these new operators, 



H = -tE -^M , (2.5) 
where the new couplings are given by 

^ = ^—(h-r(3) , r = -^{f3-sh) . (2.6) 
1 — rs I — rs 



We note that the general ansatz given by Eqs. |2.4| to p.6| does allow for the possibility 
that the new couplings r and C, do not define orthogonal directions in the space of 
the original couplings (3 and h, i.e. they need not result from a rotation of the 
couplings P and h. 

In the presence of a non-vanishing external field the line of first order phase 
transitions, Pc{h), singles out a direction which corresponds to the low temperature, 
symmetry broken part of a temperature-driven transition. The new temperature- 
like direction r thus is identified with (3c{h). In the vicinity of the critical endpoint 
the slope of this line determines the mixing parameter r, 



dh / h=hc 



(2.7) 



The second mixing parameter, s, is determined by demanding that the energy-like 
fiuctuations and those of the ordering field are uncorrelated, 

{6M ■6E)=0 , (2.8) 

with SX = X — (X) for X = M and E. This insures that the expectation value of 
the new ordering field operator, (M), fulfills a basic property of an order parameter, 
i. e. for ^ = it stays r-independent in the symmetric phase. Eqs. and 2]^ give 



two independent conditions which are sufficient to determine the mixing parameters 
r and s. 

All our simulations have been performed on lattices of size with L = 40, 50, 60 
and 70. We use periodic boundary conditions and perform simulations with a cluster 
algorithm described in Ref. |T0|. Independent of the lattice size we call a new 



configuration the spin configuration obtained after 1000 cluster updates. Typically 
we have for each value of the external field h performed simulations at 3 to 4 /3- values 
in the vicinity of the critical point. For each pair of couplings {f3, h) we generated 
about 10000 configurations which then have been used in a Ferrenberg-Swendsen 
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reweighting analysis [|r^ to calculate observables at intermediate parameter values. 
Autocorrelation times have been estimated by monitoring the time evolution of 
the energy E. Close to the pseudo-critical points they vary between 7 and 25 
configurations on the smallest and largest lattices, respectively. All errors quoted 
below have been obtained from a jackknife analysis. 



3 Determination of the critical endpoint 

The basic observables for the determination of the pseudo-critical couplings on finite 
lattices are susceptibilities constructed from the Hamiltonian H, and the magneti- 
zation M, 

CL = ;^((^')-(^)') ' (3-1) 

XL = ;^((M2)-(M)2) . (3.2) 

The location of the maxima in these observables define pseudo-critical couplings, 
Pc,L{h), which may differ for different observables. We generally find that the loca- 
tions of CL,max and XL,max differ by about one standard deviation of the statistical 
errors on Pc,L{h) for our smaller lattices and agree within errors for the largest lattice, 
L = 70. The maxima have been obtained from a combined Ferrenberg-Swendsen 
analysis which takes into account all data sets at fixed values of h at /5-values in the 
vicinity of the pseudo-critical points. In Tab. H we quote results for Pc,L{h) obtained 
from XL,max for couplings, h, in the vicinity of the critical endpoint. 

A first indication for the location of the critical region is obtained from an analysis 
of the volume dependence of the peak heights in the susceptibilities. In the region 
of first order phase transitions {h < he) XL,max is expected to increase proportional 
to the volume, while for h > he the peak heights will approach a finite value in the 
infinite volume limit. Results for XL,max obtained on different size lattices are shown 
in Fig. 1. From this figure as well as from a similar analysis of CL,max it is clear 
that the critical endpoint will be located in the interval 0.0005 < h < 0.001. In this 
interval we find that the dependence of the pseudo-critical couplings /3c,l(^) on the 
external field h is well approximated by a leading order Taylor-expansion, 

PcAhi) - PcAh2) = —ih- h^) , /ii, h2 G [0.0005, 0.001] . (3.3) 

^ L 

For our largest lattice, L = 70, this dependence is shown in Fig. 2. A straight 
line fit to these data yields, ryo = —0.689(8). Results from other lattice sizes are 
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Figure 1: The maxima of xl (a) and xl in units of the volume F = (b) as a 
function of the external field h on various size lattices. 
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{h) 




h 


L = 40 


L = 50 


L = m 


L = 70 


0.00050 


0.549800(16) 


0.549798(10) 


0.549792 (8) 


0.549790 (6) 


0.00055 


0.549742(16) 


0.549724(12) 


0.549718 (8) 


0.549716 (6) 


0.00060 


0.549668(18) 


0.549652(12) 


0.549644 (8) 


0.549642 (6) 


0.00065 


0.549596(20) 


0.549578(14) 


0.549572 (8) 


0.549568 (6) 


0.00070 


0.549524(20) 


0.549506(14) 


0.549498 (8) 


0.549496 (6) 


0.00075 


0.549454(20) 


0.549434(14) 


0.549426 (8) 


0.549424 (6) 


0.00080 


0.549382(22) 


0.549362(16) 


0.549354(10) 


0.549352 (6) 


0.00085 


0.549312(22) 


0.549292(16) 


0.549282(10) 


0.549280 (8) 


0.00090 


0.549240(22) 


0.549220(18) 


0.549212(10) 


0.549208(10) 


0.00095 


0.549170(24) 


0.549150(18) 


0.549140(10) 


0.549136(10) 



Table 1: Pseudo-critical couplings determined from the location of the peak in xl 
on different size lattices at various values of the external field h. 



L 


40 


50 


60 


70 


oo 


tl 


-0.706(21) 


-0.694(15) 


-0.690 (9) 


-0.689 (8) 


-0.685(10) 


sl 


0.696 (1) 


0.696 (1) 


0.694 (1) 


0.690 (3) 


0.690(2) 



Table 2: Slope parameter, r^, determined from the dependence of the pseudo-critical 
couplings on the external field h (Eq. p.7|) and the second mixing parameter sl 
determined from Eq. for different lattice sizes. 

summarized in Tab. |^. The slope parameter is slightly volume dependent which, 
of course, reflects the volume dependence of the pseudo-critical couplings. We thus 
have extrapolated to the infinite volume limit using the ansatz = + c/L^. 
This yields Too = —0.685(10) and fixes the mixing parameter r defined in Eq. |2.4| . 

The mixing parameter determined above may be used in connection with 
Eq. |2]^ to determine the second mixing parameter sl- We have done so and within 
statistical errors we found sl = —rL, i-e. the couplings (t, ^) can be obtained from 
a rotation in the space of the original couplings {(3,h). An alternative approach 
thus is to assume = —r^ and use Eq. |2.^ to determine s^. In this way we also 
can make use of the information on diagonal correlations ((5X)^), X = E, M and 
determine sl from a diagonalization of the fluctuation matrix 0, |^. The results 
obtained in this way for sl on different size lattice are also given in Tab. |^. We note 
that within errors this approach is consistent with the determination of from the 
slope of the critical line. The statistical errors are, however, significantly smaller. 
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0.5498 4: 



0.5496 



0.5494 



0.5492 




0.0005 0.0006 0.0007 0.0008 0.0009 



Figure 2: The dependence of the critical couphng on the external field on a 70^ 
lattice in the vicinity of the second order endpoint whose location is indicated by 
an asterisk. 



This shows that the transformation of variables, {l3,h) — > (r, ^) defined in Eq. 
indeed is a rotation. In the following we will use for the mixing parameters r = —s 
with s = 0.690. 

Having fixed the mixing parameters we can perform an analysis of the critical 
behaviour in the vicinity of the critical endpoint using standard techniques to study 
temperature driven second order phase transitions in the absence of an external 
symmetry breaking field. In particular, we can use Binder cumulants to locate the 
critical endpoint, 



B 



3,L 



((5M)2)3/2 



B. 



4,L 



(3.4) 



For given values of h and L we find that the cumulant B3^L{(3,h) vanishes and 
B4^L{l3,h) acquires a minimum at values of the coupling (3 which agree with each 
other within statistical errors. This, in turn, defines a pseudo-critical coupling, 
which again agrees within errors with the pseudo-critical couplings given in Tab. 
The minima of the second Binder cumulant, B^ l, calculated on different size lattices 
should have a unique crossing point when plotted, for instance, versus the external 
field h. This is indeed the case, as can be seen in Fig. ||. The crossing point yields 
the critical field he at the second order endpoint. 
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Figure 3: The Binder cumulant B^ l versus h for different lattice sizes. The horizon- 
tal line shows the universal value of the Binder cumulant for the universality class 
of the three dimensional Ising model, B4 = 1.604(1) [ |T^ . 

From the crossing points of the Binder cumulants we find for the critical endpoint 



/?c = 0.54938(2) , = 0.000775(10) . (3.5) 

or equivalently in terms of the rotated couplings, 

(r^,^J = (0.37182(2), 0.25733(2)) . (3.6) 

At this point the Binder cumulant takes on the value B4 = 1.609(4)(10), where the 
first error only takes into account the fluctuation of B^^l on the different size lattices 
at {jSc, he) and the second error estimates the uncertainties arising from the errors 
on the location of the endpoint. Our result agrees within errors with the value found 
for the three dimensional Ising model, S4 = 1.604(1) ||12|| . 



4 Universality class of the endpoint 



The value of the Binder cumulant B^ determined at the critical endpoint strongly 
suggests that the critical endpoint belongs to the universality class of the 3-d Ising 
model. Rather impressive support for Ising-like behaviour at the critical endpoint, 
which at the same time also demonstrates the importance of the correct choice of 
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Figure 4: Correlations between the original fields E and M (left) and the new 
energy-like and ordering field like operators E and M. Shown are contour plots 
obtained from joint histograms for these observables on a 70^ lattice at {(3c, K). The 
histograms have been obtained from a Ferrenberg-Swendsen reweighting in (3 and 
h. Actually shown are the normalized observables which have vanishing expectation 
value and a variance that is equal to unity, (X - (X))/((X - {X)fy/^. 



energy- and ordering field like variables, is given by the contour plots obtained from 
the joined probability distributions (two-dimensional histograms) of the operators 
{E,M) and {E,M), respectively. These are shown in Fig. ^ for the 70^ lattice 
at {f3c,hc)- Although the contour plot for {E,M) is still slightly asymmetric, a 
comparison with the corresponding contour plot of the 3-d Ising model ^ clearly 
shows the same characteristic fingerprint. We have checked that this asymmetry is 
well within the uncertainties of our determinations of (/3c, he) as well as the mixing 
parameters r and s. 

Further support for the Ising universality class of the critical endpoint comes from 
a more conventional finite size scaling analysis performed for the newly defined mixed 
operators M, E and the related susceptibilities. In particular we have considered 
the order parameter 



m(r,0 = ;^(M(r,0-(M(^c,ec))) 
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(4.1) 





this work 


3-d Ising 




1.60(2) 


1.587(2) 


P/u 


0.517(3) 


0.5185(15) 


7/z/ 


1.93(1) 


1.9630(30) 


{l-a)/v 


1.30(2) 


1.413(2) 



Table 3: Critical exponents determined from a finite size scaling analysis in the 
vicinity of the critical endpoint of the Potts model. The last column is based on the 



exponents for the 3-d Ising model given in [12 



and the susceptibility 



XL = Ln{m')~{\m\) 



(4.2) 



These observables indeed show the usual finite size behaviour in the vicinity of 
the critical point. For fixed ^ = C,c the susceptibility xl reaches a maximum at a 
pseudo-critical coupling Tpc. Critical exponents can then be determined from the 
scaling of the order parameter, (|m|) ~ L~^/^ , the peak height of the susceptibility, 
XL,max ~ L~^l^ ^ and the pseudo-critical couplings, (r^c — r^) ~ L^^/" . In our analysis 
of the critical behaviour we have first fixed the critical point (tc, .^c) and determined 
critical exponents from two-parameter fits. Errors have then be determined by also 
varying (tc, ^c) within the errors given in Eq. |3.6| . In Tab. ^ we summarize the results 
of this analysis and compare the calculated exponents with known values for the 3-d 
Ising model. As can be seen the agreement is quite satisfactory although in our 
analysis the lattices have not yet been large enough to reach an accuracy similar to 
that obtained for the 3-d Ising model 



Most sensitive to the correct choice of the new, mixed observables are the new 
energy-like observable, (e) = {E/L^), and the corresponding susceptibility, Ce ~ 
(de/dr)^. Their critical behaviour is controlled by the thermal critical exponent 
(yt) which fixes the critical exponent a. In particular, the energy density calculated 
at (tc,^c) is expected to scale like (e) = cq + ciL~^^~"^^'^ at (tc^c)- It is obvious 
from Eq. |2.4| that this scaling behaviour can hold only if in the mixed observable e 
the leading singular behaviour of E/L^ and M/L^, which in general is proportional 
to L^^/'^ ^ gets canceled. For the 3-d Ising model the exponent a is small (a = 0.11). 
As a consequence the dominant singular behaviour in e is expected to behave like 
~ YsXhei than ~ as would be the case for E/L^ and MjL^ separately. 

As can be seen from Tab. |^ this is supported by our analysis of the finite size 
scaling behaviour of the energy-like observable (e), although in this case the critical 
exponent shows the largest deviation from the corresponding value of the 3-d Ising 
model. 
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5 Conclusions 



We have determined the critical endpoint of the 3-d, 3-state Potts model in the 
presence of an external ordering field and have analyzed the critical behaviour in 
the vicinity of this second order phase transition. New energy-like and ordering 
field like couplings and observables have been obtained through a rotation of the 
original couplings and fields. At the critical endpoint of the line of first order phase 
transitions the joined probability distribution of these new observables shows a cor- 
relation pattern which is characteristic for the 3-d Ising model. A finite size scaling 
analysis performed with the newly defined rotated observables further supports that 
the critical endpoint belongs to the universahty class of the 3-d Ising model. 
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